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I. INTRODUCTION 

Self organization is one of the most extraordinary tools 
in Nature which assembles its elementary building blocks, 
atoms and molecules, in the form of solid substances. 
Crystals, in this sense, represent the class of a vast num¬ 
ber of solid materials in which atoms (molecules) are 
brought together in a perfect, spatially-periodic manner. 
A macroscopic crystal, is therefore made by repeating 
its smallest microscopic period - the crystal unit cell - 
infinitely in space (Fig. [^a)). Therefore, when trans¬ 
lating the crystal by an integer number of unit cell size, 
the crystal coincides with itself. This property is in the 
basis of modern solid state physics and crystallography 
and is known as the translational symmetry. Crystals are 
thus ordered materials in which atoms show both short- 
(within a couple of nearest neighbors) and long-range or¬ 
der. Due to the periodic lattice potential, the energy 
states within a crystal are extended, and the electrons 
can freely diffuse through it. 

At the opposite extreme of crystalline order there are 
the disordered or amorphous solid materials, in which 
the short-range order may still preserve, while the long- 
range order misses completely (Fig. [^c)). In these kind 


of solids, e.g. in heavily doped semiconductors or glasses, 
periodicity lacks and the atomic potential varies in a ran¬ 
dom manner. In certain cases, when the degree of ran¬ 
domness is large enough, electron wavefunctions localize 
exponentially near individual atomic cites and the diffu¬ 
sion of charge carriers vanishes (Anderson localization) 

ID. 

With the advances in X-ray diffraction tools in the 
beginning of the XX-th century, scientists learned to 
image the signatures of long-range order of crystals 
by mapping the symmetries of their reciprocal space 
lattices. The crystal symmetries for all possible arrange¬ 
ments of perfect periodicity are classified in 230 space 
groups. On contrary, the diffraction form a disordered 
solid shows cloud-like concentric rings indicating to 
a diffraction from totally random atomic lattice and, 
consequently, a lack of any symmetry or, equivalently, 
any preferential direction in the solid. 


II. IN BETWEEN PERFECT PERIODICITY AND 
COMPLETE RANDOMNESS 


In 1982 scientists found that certain “crystalline” solids 
were showing point-like X-ray diffraction patterns with 
symmetries, e.g. 5-fold and 10-fold, that did not fit any 
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FIG. 1. Schematic views of a two-dimensional (a) periodic, 
(b) a quasiperiodic and (c) a random lattice. Below each panel 
the corresponding situation for the atomic order is shown. 




















































2 


of the known 230 space group^ [42] . These solids were 
named quasicrystals and represent a large class of ma¬ 
terials in between crystalline and disordered solids. A 
quasicrystal shows a perfect long-range order but does 
not possess a periodic lattice (Fig. Bb)). Equivalently, 
the quasicrystalline lattice can densely fill the space, how¬ 
ever, it does not posses translational symmetry. 

From a mathematical point of view, two- (2D) and 
three-dimensional (3D) quasicrystals are formed follow¬ 
ing aperiodic tiling rules. Discovered by mathematicians 
in early 60’^ the aperiodic tilings can fill densely the 
plane by a combination of scaling, rotating and repeat¬ 
ing procedures of compound patterns. Similar rules allow 
to compound forms for filling the space m- 

So far, on the way towards our understanding of the 
surrounding physical world and natural phenomena, sci¬ 
entists have adopted mathematical approaches which 
point to a reduction of an existing problem into a one¬ 
dimensional (ID) one. This is because in ID mathemat¬ 
ical models can be described analytically and solved ex¬ 
actly. Higher-dimensional problems, in many cases, are 
then constructed based on a generalization of ID models 
and their solutions, while in other cases, numerical meth¬ 
ods and laborious computer simulations should run. 

In the next section, we will address principal aspects 
of ID quasiperiodicity with a particular focus on ID Fi¬ 
bonacci chains. Further, the rest of the chapter will be 
dedicated to the electromagnetic counterpart of ID Fi¬ 
bonacci structures as a relatively simplest case of the 
large class of photonic quasicrystals. 


III. ID QUASIPERIODICITY: FIBONACCI 
CHAIN 

Perhaps, the most well known and studied example of 
an aperiodic sequence of numbers is the Fibonacci se¬ 
quence [45]. The first two numbers, called the seed, are 
two I’s and the following numbers 

2, 3, 5, 8, 13, 21, ... (1) 

are obtained by summing each time the last two of the 
lisl|^ Mathematically, the recurrence rule defines the n- 
th Fibonacci number as 

Fn = Fn-i + F^-2 (2) 


^ According to the crystallographic restriction theorem the crys¬ 
tals are limited to 2-fold, 3-fold, 4-fold, and 6-fold rotational 
symmetries. 

^ Girih tilings were used in far 1200’s for decoration of buildings 
in medieval Islamic architecture. 

^ In an alternative definition, the seed numbers are 0 and 1, but 
the recurrence rule is the same. 




FIG. 2. (a) The Fibonacci numbers and the golden ratio, (b) 
Tiling of squares with sides increasing as Fibonacci numbers 
and the golden spiral, (c) The spiral pattern of Vogel repro¬ 
ducing sunflower florets, (d) The two-component Fibonacci 
words. 


with the seeds Fq = 1 and Fi = 1. The successive Fi¬ 
bonacci numbers, as shown in Fig. [^a), are intimately 
connected through the Golden ratio 

cp = lim ^ ~ 1.61803 ... (3) 

n^oo r 2 

The deterministic nature, that follows from recurrence 
rules used to create aperiodic sequences, and, in partic¬ 
ular, that of Fibonacci type, allows to construct a vari¬ 
ety of many different “forms” - strings or patterns. A 
nice visual example is the 2D tiling of squares, which are 
scaled such that their sides equal to Fibonacci numbers 
in length, as shown in Fig. [^b). 

Scientists have observed that in many plants, such 
as sunflower and cactus, the florets and stickers are or¬ 
ganized on composite spiral-like patterns following Fi¬ 
bonacci rules [Hill. A model to reproduce the sun¬ 
flower floret pattern has been introduced by H. Vogel 
M- It constitutes in generating a scatter plot in po¬ 
lar coordinates (r = a^n, 6> = 27rn/(/:?^), where a is a 
scaling factor, n is the index of the floret and ip is 
the Golden ratio (Fig. |^c)). In this model, the angle 
Og = Iip‘^ = 137.508° is the so called Golden angle, at 
which in a circle of circumference a b the ratio a/5 of 
arcs a at 27r — Oq and h at Oq equals to (a + 5)/a. 

In ID, by using recurrence rules, one may construct 
aperiodic sequences starting from seed strings. The most 
famous example is forming so-called “aperiodic words” 
by starting with seeds A and B and recurrently applying 
some rule for the next word-generations. In particular, 
in the case of a Fibonacci sequenc^ an inflation rule is 


Another intensively studied aperiodic sequence is the Thue- 
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String 

Length 

counts (A) 

counts (R) 

A 

1 

1 

0 

B 

1 

0 

1 

AB 

2 

1 

1 

ABA 

3 

2 

1 

ABAAB 

5 

3 

2 

ABAABABA 

8 

5 

3 

ABAABABAABAAB 

13 

8 

5 


TABLE I: The Fibonacci strings and relations between its 
length and counts of components. Length/counts(A)^ (f 
and counts (A)/counts(B)^ (f for long Strings. 


used 


AB (4) 

(5) 

The first six Fibonacci strings, therefore, are A, 5, AB^ 
ABA, ABAAB and ABAABABA (Fig. |^d)). As it 
is seen, this inflation rule is equivalent to adding the 
(n — 2 )th sequence to the end of the (n — l)th in order to 
generate the nth series, Fn = Fn-iFn- 2 - One particular¬ 
ity of the Fibonacci string is that, while A A is a frequent 
instance along an Fn for large n’s, other instances such 
as BB or AAA may never appear. As it will be discussed 
in the next section, this particularity appears to have sig¬ 
nificance in the context of a ID photonic Fibonacci-type 
quasicrystal. 


IV. PHOTONS IN A ID OPTICAL POTENTIAL 

The energy spectrum of a quantum particle in a semi¬ 
conductor crystal is described by extended Bloch modes 
and consists of allowed and forbidden bands. Electronic 
crystals have their analogue in the world of electromag¬ 
netic waves (photons) called photonic crystals [27]. Pho¬ 
tonic crystals are complex materials, in which a periodic 
variation of the material dielectric constant, e(r), on a 
length scale comparable with the wavelength leads to the 
formation of energy bands where the propagation of pho¬ 
tons is allowed or forbidden. 

Likewise solid materials, photonic crystals can be peri¬ 
odic, quasiperiodic and random, reflecting the variation 
of e(r) in space. This variation, in turn, can be either 
continuous (Eig. [^, top) or abrupt with the space co¬ 
ordinate (Eig. [^a), bottom). In the first case a smooth 
variation of e(r) between ei and 62 may reflect a con¬ 
tinuous change in material composition (an Ala,Gai_a,As 
alloy, for example). In the second case an abrupt change 
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FIG. 3. (a) The variation of the dielectric constant in space 
on a length scale comparable to the light wavelength lays 
in the basis of photonic crystals. e(r) may vary either in a 
smooth or a step like manner, (b) With a step-like change of 
e(r), by using to different materials A and B, results in the 
formation of a periodic (top), a quasiperiodic (middle) or a 
completely random (bottom) optical potential for photons, 
(c) The Fourier transformation spectra of the correspond¬ 
ing bi-component layer sequences show a single sampling fre¬ 
quency (periodic), a self-similar frequency distribution pat¬ 
tern (quasiperiodic) and an uncorrelated flat spectrum with 
no any specific frequency (random). 


in e(r) is achieved when two different materials with di¬ 
electric constants ei and 62 are stack together. 

Throughout this chapter, we will focus on the optical 
properties of ID photonic quasicrystals of Fibonacci type 
which are formed by stacking two different materials, say 
A and 5, following the inflation rule discussed above. 
Prior to this, we will introduce the basic properties 
of ID photonic multilayered structures on the simple 
examples of a Distributed Bragg Reflector (DBR) and a 
Fabry-Perot microcavity. 


Distributed Bragg Reflector 

Otherwise known as a Dielectric Mirror, the function¬ 
ality of a DBR is based on the interference between elec¬ 
tromagnetic waves (EMW) which are multiply reflected 
and refracted at each boundary of the multilayer struc¬ 
ture. An EMW of frequency ujq is efficiently reflectec^ 
from a DBR if the multilayer is constructed by alternat¬ 
ing quarter-wave layers of two different materials A and 
B with refractive indices ua = \f^A ~ 

respectively. The quarter-wave condition implies that 


Morse one which uses an {A —)► AB, B -A BA} inflation | 51| . 


^ For simplicity, we will consider the case of a normal incidence in 
the following. 
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FIG. 4. The multilayer structure aud the trausmissiou spectra 
of (a) a Distributed Bragg Reflector aud (b) a Fabry-Perot 
microcavity. 


= TiBdB = Ao/4, where (Ia and (Ib are the layer 
thicknesses, Aq = 27tc/ujq is the wavelength and c is the 
EMW velocity in free space. The portion of the EMW 
intensity reflected from the multilayer embedded in air 
from both sides is given by 


R = 


^2N 


-1 2 


^2N 

J 


,2N 


^2N 

I'^A 


(6) 


where N is the number of periods {AB pairs) of the DBR. 
In absence of losses (scattering, absorption etc.) the 
transmission of the mirror is T = 1 — R. As shown in the 
example of Eig. |^a), the transmission vanishes within a 
bandwidth of Auj around the frequency uq, which means 
that in this frequency window - the photonic band gap 
- EMWs cannot propagate through the structure. The 
relative width of the photonic band gap is given approx¬ 
imately by 


Auj 4 

- = — arcsm 

CJo TT 


f \nA - 

\nA+nB ) 


( 7 ) 


Erom Eq. ^ we see that the reflectivity of the di¬ 
electric mirror grows with either increasing the number 
of periods N or increasing the refractive index contrast 
‘ crucial for making a wider pho¬ 

tonic band gap, as follows from Eq. 0 . while increasing 
N affects only the steepness of the band edges. 

A common mathematical approach to calculate the 
spectral properties of ID multilayered photonic struc¬ 
tures is based on the transfer-matrix method [37]. This 
method consists in relating the electromagnetic fields on 
two boundaries (left and right, hereafter) of the dielectric 
layer by means of the transfer matrix, T. In the case of 
the electric filed E this relation is 


Ei = M X Er 


cosJ i-sinS 

n 

i-sind cos 

c 


Er 


(8) 


where n, d and 6 = are the refractive index, the 

thickness of the layer and the phase change across it, 
respectively, Ei^r are the electric fields at the left and 


right boundaries of the layer. Equation 0 in the same 
form holds also for the magnetic field. 

Erom continuity requirements it follows that the trans¬ 
fer matrix M for a system of N layers is obtained by 
simple matrix multiplication of M’s of single layers, 
M = n with Mj being the transfer matrix of the 
jf'th layer. The complex reflectance, r{uj) and transmit¬ 
tance, t{uj) coefficients of the multilayer can be thus cal¬ 
culated from the total transfer matrix M by assuming a 
unity (normalized to the input) intensity impinging on 
the structure from one side (left) and no input from the 
other (right): 

riuj) = -, (9) 

limii + 7z7rmi2 + m2i + 7r^22 


t{uj) = -, (10) 

7 /mii + ^armi2 + rn2i + 7r^22 

where 7 / and 7 ^ are the inverse of light velocities in the 
input and output media, respectively. 

The complex quantities contain important information 
about the photonic properties of a multilayer. Eor exam¬ 
ple, the power reflection and transmission coefficients are 
calculated as i? = \r{uj)\'^ and T = |t(u;)p, respectively. 
Moreover, the amplitude and the phase of the transmit¬ 
ted light can be retrieved from t{uj). The phase 0, in 
particular, which is the argument of the complex trans¬ 
mittance, arg[t(cj)], allows one to calculate both the den¬ 
sity of states (DOS) of the photonic structure, p{uj), and 
the group velocity of propagating light, Vg{uj), as 

p{uj) = Vg{U})~'^ = ( 11 ) 

with Ltot being the total physical thickness of the mul¬ 
tilayer. Eigure |^a) plots the DOS of the same dielec¬ 
tric mirror of Eig. [^a). It shows, in particular, that the 
DOS, defined typically as the number of states (photonic 
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FIG. 5. The density of states (DOS) and the group velocity 
for the cases of (a) a dielectric mirror DBR and (b) a Fabry- 
Perot microcavity. In the last case a strong enhancement 
of the DOS and signihcantly reduced group velocity can be 
observed at the resonant frequency uj/ujq = 1. 

























































5 


modes in this case) per unit energetic (frequency) inter¬ 
val, is largely enhanced at certain frequencies adjacent to 
photonic band gap, known as band edge states. 

Following the definition of Eq. 0, the group velocity 
of a wavepacket centered at the band edge frequencies 
and propagating through the multilayer is strongly 
reduced for those spectral components which lay within 
the narrow high-DOS lines (see Fig. |^b)). As a 
consequence, the temporal evolution of a transmitted 
through DBR light pulse shows important wavepacket 
broadening as a result of a pronounced delay of these 
spectral components (mode localization). 



FIG. 6. A photonic quasiperiodic potential for electromag¬ 
netic waves is formed by alternating according the Fibonacci 
inflation rule two different dielectric layers of appropriately 
chosen optical thicknesses Uidi. 


Fabry-Perot microcavity 

The Fabry-Perot (FP) microcavity is a dielectric mir¬ 
ror with a defect. The simplest way to introduce a de¬ 
fect in a dielectric multilayer is to eliminate one of the 
Ao/4 layers of either type A or B materials. By doing 
this, two quarter-wave layers of the same material ap¬ 
pear next to each other, forming thus a Ao/2 layer which 
leads to a perfect constructive interference of transmit¬ 
ted at ujo EMW. As a result, a sharp (unity) transmission 
channel is formed at the center of the photonic band gap 
(Fig. I^b)). The narrowness of the transmission line de¬ 
pends on the reflectivity (i.e., the strength) of the dielec¬ 
tric mirrors sandwiching the half-wavelength defect layer. 
As for the reflectance of the FP microcavity, it vanishes 
exactly at the resonant frequency, while for the rest of 
frequencies shows all the spectral features of a DBR. 

It is important to note the dramatic increase of the 
DOS (Fig. [^c)) and, consequently, the suppression of 
group velocity at the microcavity resonance frequency 
(Fig. I^d)). These are signatures of a strongly localized 
defect state for which the mode wavefunction decays ex¬ 
ponentially outside the defect. The ability to confine 
EMWs in a tiny volume such as that of a defect layer, 
makes microresonator devices excellent tools for studying 
light-matter interactions at a nanometric scale. Lasing 
and Purcell enhancement of the radiative rates of cavity- 
embedded emitters are typical phenomena which largely 
exploit the confining properties of microresonators. 

Finally, after the introduction to the physics of a di¬ 
electric mirror and a microcavity, we now have all the 
basic ingredients necessary to start with the optical prop¬ 
erties of a ID photonic quasicrystal. In the last part of 
this chapter we will see how a ID Fibonacci quasicrys¬ 
tal is formed and in which way the quasiperiodic order 
modifies the optical properties of an otherwise ordered 
photonic structure. 


V. PHOTONS IN A ID QUASIPERIODIC 
POTENTIAL 

As briefly cited in the beginning of the previous sec¬ 
tion, an optical potential for electromagnetic waves can 
be realized by properly stacking thin quarter-wave layers 
made of two different dielectric materials A and B. In 
the two extremes, a strictly periodic alternation of these 
layers will form a perfectly ordered ID photonic crystal 
(DBR), while a complete randomization of the stacking 
order will give rise to a disordered photonic structure 
(Fig-ib)). 

In terms of the positional order of A and B layers, the a 
ID quasicrystal stands in between the perfect DBR and a 
random photonic structure. As already introduced in the 
Section [j following the aperiodic string generation rule, 
we now can construct the photonic quasicrystal struc¬ 
ture of Fibonacci type (FQ, hereafter) in which the con¬ 
stituent A/4 layers of A and 5-type materials are stack 
together obeying F/ = F^_iF /_2 for the nth Fibonacci- 
generation [HUS]. Thus, the number of A/4 layers of 
an nth-generation FQ will correspond to the nth-order 
Fibonacci number (see Fig. [^. 

The differences between the periodic, aperiodic and 
completely random arrangements of two-component ID 
systems of the same length (all A’s and 5’s together) can 
be recognized by looking at the Fourier transform spec¬ 
tra of these sequences. Figure c) reports such spectra 
for a DBR, a 15th order FQ and a positional-random 
multilayer stack, all containing 987 layers. We see imme¬ 
diately that in the case of the DBR (top panel) a single 
occurrence frequency is present, indicating the perfect al¬ 
ternation of A and B layers. Instead, the flat spectrum 
of the random structure (bottom panel) is a strong evi¬ 
dence that no any correlation is present in the sequence 
of alternation of layer types. More interesting and rich 
is the Fourier transformation spectrum of the FQ (mid¬ 
dle panel). We notice that this A/5-appearance spec¬ 
trum represents a self-similar structure full of intense 
and weaker peaks - signatures of long-range order and 
absence of perfect short-range periodicity. 
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Electronic energy spectrum of ID FQ 

Various authors have studied the energy spectrum and 
the properties of ID quasicrystals. Early studies were 
mainly focused on electronic quasicrystals and much 
work has been done there. In particular, it has been 
reported that the typical structure of the energy spec¬ 
trum, T{uj), of a Fibonacci system has a self-similar na¬ 
ture [28l [29] . More technically, it is a Cantor set with 
zero Lebesgue measure [28|. In other words, an arbi¬ 
trary chosen frequency lays within the gap with a unity 
probability, and the numerous gaps are densely filling the 
spectrum. In these spectral gaps, known as pseudo band 
gaps^ the DOS tends to zero P [34]. 

It was further shown that as a consequence of the 
multi-fractal nature of the energy spectrum, the wave- 
functions of a Fibonacci quasicrystal show certain exotic 
localization properties [29]. Namely, the wavefunctions 
are neither localized nor extended - is is said they are 
critically localized - and thus decay at a rate slower than 
an exponential one nanHiiiH]. These statements have 
been multiply confirmed through numerical calculations 
and experiments [26| [31] |39[ [47[ [50] [55] . Starting from 
mid 90’s different studies with their optical (photonic) 
counterparts confirmed major analogies to the electronic 
case HEl. 

Energy spectrum of ID photonic FQ 

From electronic to photonic - Interference is a wave 
phenomenon [43]. The wave nature of both electrons 
and photons is the basis of important analogies between 
electronic charge transport phenomena in solids and light 
propagation in complex dielectric systems [46] . The spec¬ 
trum of these analogies includes but is not limited to phe¬ 
nomena such as optical counterpart of weak localization 
[301E31E7], Anderson localization [SlIlIlISS], short and 
long range correlations m E2, and universal conduc¬ 
tance fluctuations lu, photonic Bloch oscillations m 
and Zener tunneling of light waves [23] . 

Since photons are uncharged particles and interact 
extremely weakly with each other at relatively low light 
intensities, a light wave packet, that propagates through 
a dielectric system, remains coherent for much longer 
times than charged particles. This means that dynamic 
interference effects could be isolated and studied more 
easily with light than with electrons. In view of this, ID 
photonic Fibonacci quasicrystals were not only realized 
and spectrally characterized in static experiments, 
but more importantly, photon wavepacket propagation 
through optical states around the pseudo band gap of 
the structure was studied thoroughly Elia]. 

Transmission spectrum - In general, the energy spec¬ 
trum, E{uj)^ of a photonic structure is often identified 



pseudo-gaps 




FIG. 7. (a) The transmission spectrum of a 13th order FQ 
with 377 layers. The inset shows the zoom around the central 
frequency revealing the fine modal features, (b) The evolution 
of the transmission spectra with increasing Fibonacci order: 
the pronounced deepening of low-T regions is accompanied 
by the appearance of a richer resonant structure, (c) The 
shadowed regions in the spectrum of an F9 structure indicate 
the so-called pseudo band gaps, (d) The same is illustrated 
on the DOS spectrum. 

to its transmission spectrum T{uj). We will use this last 
definition throughout the rest of the chapter. As men¬ 
tioned in the Section [| the transmission spectrum of a 
multilayered FQ can be easily calculated according the 
transfer-matrix formalism (Eqs. (§ and ([To])). 

Figure [^a) shows the transmission spectrum of a 
377-layer 13th-order Fibonacci sequence with chosen 
refractive indices of = 3 and = 1. This spectrum 
is different form that of a perfect DBR but it represents 
a mixture of general spectral features of a dielectric 
mirror and a FP-microcavity in terms of clearly vis¬ 
ible gap regions and narrow transmission lines. The 
calculated T(ci;)-spectra for four successive Fibonacci 
series are plotted in Fig. [^b). The logarithmic scale of 
transmission shows how the spectra are evolved from 
F 4 to Fj reflecting the deepening of prohibited spectral 
regions and the appearance of new narrow spectral 
features. 


Experiments with ID photonic FQ’s 

One of the first realizations of ID FQ’s and the spectral 
characterization of samples has been reported by Geller- 
mann et al in 1994 m- An electron-gun evaporation 
technique was used to grow up to 55 layer (F9) sam¬ 
ples composed by silica and titanium dioxide layers al¬ 
ternating following the Fibonacci rule. The quarter-wave 
(Ao/4) layers were centered a wavelength of Aq = 700 nm 
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and static transmission spectra of the various Fibonacci 
samples were recorded around Aq. The experimental re¬ 
sults were successfully compared to the theory based on 
the transfer-matrix approach. As a typical drawback 
from the control over growth parameters, in particular, 
a variation of layer refractive index with increasing the 
structure size, certain deviation from the nominal optical 
path {ui X di) was observed too. By fitting the experi¬ 
mental experimental spectrum a 5% of refractive index 
drift was estimated by the authors for F9 samples. 

In works by Dal Negro et al [8] and Ghulinyan et. 
al [21] dielectric Fibonacci quasicrystals were realized 
by electrochemical etching of silicon. Porous silicon has 
interesting optical properties [7] and offers an excellent 
platform for testing ideas experimentally in ID multilay¬ 
ered photonic structures 0 Hums SO]- When grown in 
highly doped p-type silicon, porous silicon behaves as an 
optically homogeneous dielectric material with an effec¬ 
tive refractive index n determined by its porosity. Since 
it is fabricated by a self-limiting electrochemical process, 
one can control the refractive index of each layer by vary¬ 
ing the electrochemical current during the fabrication. 
Typically, this last follows a step variation between two 
constant values, and therefore results in a step-like vari¬ 
ation of the effective refractive index profile through the 
whole structure (the bottom panel case of Fig.j^a)). The 
thickness d of an individual porous layer is then deter¬ 
mined by the etching duration under a constant current 
density. Hence, the control of etching current and the 
time guarantee a careful control over the optical path 
throughout the multilayered sample. 

This growth technique was used successfully to real¬ 
ize porous silicon-based ID FQ’s on silicon substrates 
and to characterize them both in static transmission 
and ultrashort-pulse interferometric experiments at near- 
infrared (NIR) wavelengths [8]. In particular, mode beat¬ 
ing, sizable field enhancement, strong pulse stretching, 
and a strongly reduced group velocity in the band edge 
region around a Fibonacci pseudo band gap were all ob¬ 
served and interpreted in terms of Fibonacci band-edge 
resonances. 

A step forward in this direction was preformed by fur¬ 
ther improving the electrochemical growth technique [21] , 
which allowed in particular to realize free-standing mul¬ 
tilayered FQ’s (without Si substrate). By using free¬ 
standing structures, in general, natural drifts of the nom¬ 
inal optical path can be controlled at a higher precision 
permitting to realize thick structures with high optical 
quality [20] [22]. Moreover, with such samples not only 
simpler static transmission experiments become accessi¬ 
ble but as well one avoids possible large time-offsets due 
to the silicon substrate during pulse propagation stud¬ 
ies. Along with the observation of all the various effects 
reported in the previous study [8], experiments with the 
new samples provided a further insight into the physics 
of photonic ID FQ’s. 


Experiment F9 (55 layers) 

- TMM-calculation 



Experiment F12 (233 layers) 

- TMM - calculation 
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FIG. 8. (a) The transmission spectrum of the porous silicon- 
based F9 sample is compared to the transfer-matrix calcu¬ 
lations with 1% of optical path gradient, (b) The same for 
the F12 sample. The gradient is of 4% in this case, (c) The 
ultrafast pulse transmission setup based on an optical gating 
technique, (d) Pulse propagation experiments were performed 
over the spectral range, which corresponds to the shaded area 
of panel (b). Dotted lines show examples of Gaussian profiles 
of the ultrashort pulses centered at different wavelengths m- 


In Figs. [^a),(b) the measured transmission spectra of 
F9 and FI2 (233 layers) are shown. The transfer-matrix 
calculations, considering the presence of possible drifts 
in growth parameters, revealed a 1% and 4% of optical 
path gradient (linear to first approximation) for F9 and 
FI2 samples, respectively (solid lines Figs. [^a),(b)). 
The central wavelength of the F9 sample was chosen to 
be approximately 1500 nm’s. Instead, the sample F12 
was intentionally realized such that the edge states of 
one of the pseudo band gaps are spectrally positioned 
around 1500 nm (see the shaded area in Fig.j^b)), where 
ultrashort pulse propagation studies could be performed 
(available range for sweeping the pulse wavelength). 

Ultrashort pulse propagation through ID FQ’s - Ultra¬ 
fast time-resolved transmission measurements were per¬ 
formed using an optical gating technique based on sig¬ 
nal upconversion (Fig. |^c)). While being simpler com¬ 
pared to an interferometric system [8], it is sensitive to 
the transmitted intensity only and not the phase. In 
experiments wavelength-tunable pulses between 1400 to 
1570 nm and of 220 fs duration (bandwidth of 14 nm) 
were sent through the FI2 sample probing thus the band 
edge states (Fig. [^d)). The transmitted signal was then 
mixed in a non-linear BBO crystal with a time-delayed 
reference pulse in order to generate a sum frequency sig¬ 
nal. This last, selected with a prism and detected with a 
photodiode, is proportional to the temporal overlap be¬ 
tween the signal and reference pulse. Therefore, the time 
profile of the transmitted signal was mapped by varying 
the delay of the reference (temporal resolution of 260 fs). 

The left panel of Fig. [^reports examples of the tempo- 
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FIG. 9. Left panels show examples of the time-resolved signals 
transmitted through the FI2 sample. In the top panel the 
reference (undisturbed) pulse profile is shown, whereas the 
rest of the panels plot the transmitted pulses at wavelengths 
of 1470 nm, 1492 nm, 1500 nm, and 1510nm, respectively. 
The concept of the center of the mass of transmitted pulses 
and the corresponding delay time are explained too. Right 
panels report (a) the samples transmission spectrum, (b) the 
delay time, (c) the relative group velocity derived from the 
delay time and (d) the decay time of the signals as a function 
of the probe pulse wavelength |21j . 


FIG. 10. (a) The transfer-matrix-calculated transmission 

spectrum (solid line) of the F12 sample and the corresponding 
DOS (dotted line), (b) The experimentally measured delay 
times of the transmitted pulses (scatter data) are compared 
to the calculated ones (line). An example of the time-resolved 
transmitted intensity of a probe pulse centered at 1497 nm is 
plotted in the inset, (c) and (d) show the light intensity dis¬ 
tribution inside the FI2 sample for the ideal case (no optical 
path gradients) and for a 4% gradient, respectively. The ar¬ 
rows in the lower graph indicate the wavelength above which 
the sample transmission becomes almost zero, while internal 
modal structure can be clearly identified [21] . 


ral evolution of selected four transmitted pulses centered 
at different wavelengths of the band edge. The top panel 
shows the free-space pulse (with no sample) the maxi¬ 
mum of which is taken to be the zero of the time axis. 
An analysis of the transmitted signals at different wave¬ 
lengths provides with rich information about the nature 
of the band edge states; 

Disorder-related effeets in ID FQ^s - From a quick 
look at the calculated DOS and the transmission shown 
in Fig. [lQ|(a) it is quite evident that while the first shows 
peaks until a wavelength of 1550 nm, the second appears 
to vanish already above 1520 nm. This fact indicates 
that there do exist states in the spectral range between 
1520 nm and 1550 nm, but they barely transmit the light 
for some reason. The origin of this phenomenon comes 
from the presence of a certain degree of disorder in the 
quasicrystal and will be discussed in the following. 

A comparison between transmission spectra, DOS, and 
group velocity of a perfect DBR, an ideal Fibonacci, and 
a Fibonacci structure with certain disorder has been re¬ 
ported and discussed in [35]. In an ideal non-absorbing 
structure, in which the optical path n^di is exactly the 
same for each individual layer, whenever the DOS shows 
a peak, it has a unity transmission. Consequently, both 
the DOS and T{uj) vanish in the band gap region. 

The disorder inside a ID photonic structure can be 
either uncorrelated to the position inside the structure or 
it can manifest as a continuous drift of the optical path 
(a linear gradient to first approximation). This last is 
the typical situation when during the preparation of the 


multilayered samples the growth rate (layer thickness) 
or the material composition (refractive index) is slowly 
varying through the final structure [i HUSH Eg. We 
will therefore examine the effect of this type disorder. 

When a drift is present, the low transmission regime 
apparently extends for a wider wavelength range than 
the DOS gap. The reason can be found by examining 
the calculated scattering states maps of the Fibonacci 
sample, reported in Fig. [lQ|^c),(d). The scattering states 
map for the ideal case with fiat photonic band structure 
is shown in the top panel. The linear drift in the optical 
path, d{nidi)/dx^ along the depth x in the multilayered 
structure acts as a built-in bias that tilts the band edge 
(see Fig. [Toj^d)). This band tilting behaves as the ana¬ 
logue of an external static electric field applied to an elec¬ 
tronic superlattice structure. Such controlled realization 
of optical path gradients was exploited recently to ob¬ 
serve time-resolved photonic Bloch oscillations [40] and 
Zener tunneling of light in optical superlattice structures 
with tilted photonic band structure [2g. Moreover, re¬ 
searchers have also managed to change the degree of band 
tilting by means of vapor flows through porous multilay¬ 
ers [31IT6]. 

Generally, and, in particular, in the 12th order Fi¬ 
bonacci sample, the drift-induced band tilting reduces 
the spatial extension of the optical modes which are 
closely situated at the edge of the pseudo band gap. As 
a consequence, the first band edge states do not extend 
over the whole sample and, therefore, their contribution 
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Flat bands 


Tilted bands 



FIG. 11. (Left graph) In the absence of optical path gradients 
the photonic bands are flat and the photon states are extended 
throughout the whole structure. (Right graph) The effect of 
an optical path gradient on photonic bands is in close analogy 
to the effect of a static electric field on energy bands in an 
electronic crystal. In the photonic case the band edge states 
are no more extended Bloch states and are confined spatially 
between the physical edge of the sample and the photonic 
band gap. Thus, photon transmission through these states is 
largely suppressed. 


to high transmission channels reduces significantly. In 
other words, a photon may be injected easily into such 
a spatially squeezed state from one side of the sample, 
however, in order to be transmitted it needs to tunnel 
through the band gap which extends for the rest of the 
sample (see Fig. 11). This explains the fact that the 
transmission spectra and density of modes differ in the 
region of the pseudo band gap. Equivalently, the long de¬ 
lay times observed around the edge of the pseudo band 
gap can be associated with the time needed to propagate 
through the band edge states. The long decay times are 
directly related to the spatial confinement of these states. 


Further numerical calculations in this direction have 
been also performed. In particular, the intensity of the 
electric field distribution inside the Fibonacci structure 
for few band edge states was investigated for both an 
ideal sample and one with optical path gradient. The re¬ 
sults of these calculations are shown in Fig. in which 
also the DOS at the band edge is reported in top pan¬ 
els. In the absence of optical path gradients (left panels), 
the distribution of the electric field intensity shows the 
characteristic self-similar structure of the band edge res¬ 
onances of a Fibonacci quasicrystal. 


It can be seen from the right panels of Fig. 12 that 
while the main features of the intensity distribution pre¬ 
serve when a drift is introduced (4% in this case as in 
the experiment with F12 samples), the profiles become 
asymmetric with respect to the center of the sample. 
Quite visibly, the intensity pattern squeezes more and 
more towards one side of the quasicrystal as the mode 
wavelength approaches the pseudo band gap. In accor¬ 
dance with the situation shown in Fig. [^for the tilted 
band case, the closest to the gap states are spatially more 
confined between the sample surface and the tilted band 
gap. Accordingly, when approaching the band edge (from 
panel A to panel B of the tilted band case), the inten- 
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FIG. 12. The electric field intensity distribution inside the 
multilayered sample for the first four band edge states in an 
ideal Fibonacci structure (left panels) and a 4% drifted one 
(right panels). The corresponding DOS spectra around the 
band edge are shown in the upper panels. 


sity profile distorts heavily and the electric field decays 
efficiently through the rest of the structure. Thus, these 
calculations show that the natural drifts, occurring dur¬ 
ing the growth of the ID photonic FQ, alter the spatial 
extension of the first few band edge modes, but do not 
destroy their specific structure. Higher order band edge 
modes are less affected by the optical path gradients. 


Origin of band edge states and pseudo band gaps 

In the previous subsection we reviewed the recent 
progress in photonic ID Fibonacci quasicrystals in terms 
of their realization, spectral characterization and ultra- 
short pulse propagation studies. We learnt that these 
photonic structures are somewhat in between ideally pe¬ 
riodic photonic crystals and disordered complex dielectric 
structures. When describing them, terms such as pseudo 
band gaps and band edge states with exotie (eritieal) lo- 
ealization properties are often used. Even though, the 
reader might still have it unclear why “pseudo” is used 
or what do we mean by saying “critical localization” and 
why does it happen? The scientific literature is rather 
succinct in this sense and rigorous mathematical descrip¬ 
tions are not easy to interpret. 

Through the rest of this chapter, we will offer the 
reader a new interpretation of these properties for the 
ID photonic Eibonacci quasicrystals. We will show, in 
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Quantum wells 
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FIG. 13. In analogy to the electronic coupling of separate 
quantum wells in a semiconductor superlattice (a), an optical 
superlattice can be realized when optical cavities are brought 
together (b), which results in the formation of a miniband of 
extended photonic states, (c) A scanning electron micrograph 
of a one-dimensional porous silicon optical superlattice with 
seven coupled microcavities. 



Frequency, co/coo Frequency, m/coo 


FIG. 14. Broad range transmission spectra of (a) a dielectric 
mirror, (b) a Fabry-Perot microcavity and (c) five coupled 
microcavities. The zooms of spectra around the band edge 
states are shown in relative panels in the right. 


particular, how the particularity of the inflation rule of 
Fibonacci-type affects the properties of a ID photonic 
FQ and in which way the so-called pseudo band gaps 
and the band edge states with non-exponential decay 
are formed? Perhaps, we will be rather surprised to 
realize that the Fibonacci quasicrystal is a particular 
arrangement of simple photonic building blocks - DBR’s 
and FP microcavities - arranged all together to form a 
coupled microcavity system. 

Coupled microcavity optical superlattice - In a close 
analogy to a semiconductor double-quantum well, the 
optical coupling between identical cavities lifts up the 
mode degeneracy inducing a repulsion, which splits the 
single state at a frequency ujq into two new ones at fre¬ 
quencies cjo ± Acj [2j [36] |38l [49] . Degenerate mode cou¬ 
pling between N cavities therefore leads to the forma¬ 
tion of a miniband of N new optical states, which are 
densely packed around the resonant frequency. This way, 
in a very same analogy to the electronic superlattices 
(Fig.[^a)) m, an optical superlattice can be made by 
coupling degenerate optical resonators (cavities) within 
the same photonic structure (Fig. [^b))[nii2ni[22i[a2i 

Eg. 

In one dimension, an optical superlattice can be re¬ 
alized by stacking two different quarter-wave dielectric 
layers A and 5 in a way to form identical cavities sep¬ 
arated by dielectric Bragg mirrors (Fig. p^c)). In a 
generic form, the layer sequence of such a multilayer stack 
looks like [{BA)r^B] (AA)i [{BA)r^B] {AA)2 ... (AA)^ 
[{BA)nB], where a series of m microcavities (AA)^ are 
coupled to each other through the [(^A)^^] DBRs of 
(n+1/2) periods. The amount of mode splitting of cavity 
resonances is given by the strength of the coupling mir¬ 
rors. As described in the Section [j for a given refractive 
index contrast, the reflectance of the DBR improves with 


increasing the number of periods. Similarly, for a fixed 
number of periods, the mirror reflectance grows with an 
increase in the refractive index contrast, i.e. the stronger 
the mirrors reflect, the weaker is the mode repulsion. 

Figure [T^ shows the formation of a characteristic trans¬ 
mission spectrum of an optical superlattice of five cou¬ 
pled microcavities. First, let us look at the spectrum of 
a single DBR which has a layer sequence [{BA)nB] with 
n = 5 (top panel)|^ As discussed previously it is char¬ 
acterized by a photonic band gap centered at a relative 
frequency ujq/uj = 1 and extends to a certain frequency 
range around it (indicated by the shaded area). 

Next, we introduce a defect in the DBR layer se¬ 
quence creating thus a A/2 cavity sandwiched between 
two DBR’s: this creates a sharp transmission peak - the 
Fabry-Perot resonance - around the central frequency 
and splits the DBR band gap in two forbidden spectral 
regions (Fig. p^b)). It is important to underline that 
the new band gaps on the left and right sides of the FP 
peak are part of the same and unique band gap of the 
DBR sequence [BABA .. .B] and therefore have the cen¬ 
ter frequency ujq/uj. 

Moreover, the FP-resonators’s spectrum contains other 
details which are absent in the DBR’s spectrum. While 
these details are not evident at a first glance, they can 
be revealed by looking at the blow up of the correspond¬ 
ing T{uj) curves in the vicinity of the band edge. From 
the right panels of Fig. we notice, in particular, that 
with respect to the DBR’s spectrum the band edge states 
of the FP structure have undergone spectral reshaping. 
Namely, a paired (double-peaked) band edge state is ap- 


The exact number of DBR periods n is not essential throughout 
this analysis. 
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pearing now. Understanding the origin of these doublet 
is important as it will be a crucial point when analyzing 
the origin of Fibonacci pseudo band gaps. 

Recalling the layer sequence of the ID FP microcavity: 


(a) 


4 coupled microcavities - 21 layers 


m :.. m :.: XU :. XU 



DBR1 DBR2 DBR2 DBR2 DBR1 


[BABA ...B]AA [BABA... B] 

and remembering that by definition they have the same 
central frequency, we notice that two identical DBR’s 
are separated through a layer A A, which can be inter¬ 
preted as a system of coupled DBR’s! As already dis¬ 
cussed above, such a coupling should induce degenerate 
mode splitting. Thus the origin of the band edge doublet 
comes from a coupling of the sandwiching DBR’s in the 
FP structure sequence. This coupling is valid for all DBR 
states. In fact, by looking back at broad range transmis¬ 
sion spectrum of the FP microcavity (Fig. p^b)), we 
notice that the doublets are present for all other modes 
far from from the band edge. Moreover, the formation of 
these doublets induces relatively deep low-transmission 
(gap-like) regions between neighboring duplets. The suc¬ 
cessive gap-like regions become shallower when moving 
further from the band edge (shallower shaded areas in 
Fig.l^b)). 

These observations become more evident when looking 
at the spectrum of a five coupled microcavity structure 
(Fig. [l^c)). First, we notice the formation of photonic 
miniband states, five in total, which are densely packed 
around ujq. As in the case of the single FP-structure, 
this miniband is separating the fundamental DBR band 
gap in two regions at uj < ujq and uj > ujq. Secondly, 
from corresponding zoom of the band gap edge we notice 
that now we have a multi-peaked modal structure (right 
bottom panel). Importantly, one counts six peaks which 
correspond exactly to the total number of DBRs in a five 
coupled microcavity structure. Moreover, the secondary 
gap-like features become much deeper and well defined. 
Finally, these observation confirm the fact that in anal¬ 
ogy to the coupling between identical A/2 AA layers also 
the dielectric mirrors couple within the same photonic 
crystal. 

It is now the right moment to look at the formation 
of the band structure of a ID photonic Fibonacci 
quasicrystal. 

ID photonic Fibonacci quasicrystal - Let us first ex¬ 
amine the layer sequence of the FQ. As an example we 
consider the 7th order Fibonacci sequence, which reads 

ABAABABAABAABABAABABA. 

Now, remembering that in the photonic case (i) A and B 
layers have a quarter-wave thickness and, hence, (ii) an 
A A layer forms a A/2 cavity, we can write this Fibonacci 
sequence in the form 

[AB] AA [BAB] AA [B] AA [BAB] A A [BABA]. 



Scattering States Map, i-og(|Ep) 





FBG2 



0 Transmission 


1 


FIG. 15. (a), (b) The schematics of layer alternation in a four 
coupled microcavity and in a 7th order FQ structure, (c) The 
scattering states map and (e) the corresponding transmission 
spectrum of the coupled microcavity structure, (d) and (f) 
show the same quantities for the ID FQ. 


We notice immediately, that this Fibonacci multilayer 
is nothing else than a system of four coupled mi¬ 
crocavities with unequal dielectric mirrors (see 
Fig.pl^a) and (b)). It is important to note here again the 
fact that due to the recurrence rule of the Fibonacci type, 
the A A is a frequent instance along an Fn for large n’s, 
while other instances such as BB or AAA may never ap¬ 
pear. In view of this, we can say now that a ID photonic 
Fibonacci quasicrystal, independently on its recurrence 
order n, will be a coupled A/2-microcavity system with 
unbalanced DBR’sH] 

For comparison, in Fig. [T^c),(d) the light intensity 
distributions inside a four coupled microcavity (4CMC) 
structure with balanced DBR’s and an F7 ID pho¬ 
tonic quasicrystal are shown. Their respective transmis¬ 
sion spectra are plotted in Fig. [T^e),(f). We note that 


^ In order to form coupled cavities (three at least, since two cannot 
be formed), n should be larger than 6. F4 and F5 have only one 
cavity, while no cavity is formed in lower orders FI, F2, and F3. 
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ID FQ 

counts(A) 

counts (R) 

counts(AA) 

ABA 

2 

1 

0 

AB\AA\B 

3 

2 

1 

AB\AA\BABA 

5 

3 

1 

AB\AA\BAB\AA\B\AA\B 

8 

5 

3 

AB\AA\BAB\AA\B\AA\BAB\AA\BABA 

13 

8 

4 


TABLE II: The occurrence of AA-cavities in layer sequences of ID FQs (from 3rd to 7th Fibonacci orders). 


the weaker external mirrors in the case of the 4CMC 
structure are chosen intentionally such to count a total 
of 21 layers as in a F7 structure. It is however important 
that the intra-cavity (coupling) mirrors are all identical. 
These intensity maps visualize nicely the layer positions 
of the A/2-cavities. In particular, one can observe that in 
the case of 4CMC structure they are evenly distributed 
through the multilayer (see the bright intensity spots 
around cj/cjo = 1 ) in contrast to their inhomogeneous 
distribution in an F7 quasicrystal. 

It follows thus that the optical coupling of cavities in¬ 
side a Fibonacci multilayer is depth-dependent. Also, it 
is clear that this coupling is not changing in a random 
manner with the position inside the multilayer, but is re¬ 
lated to the exact DBR sequences which strictly reflect 
the Fibonacci string generation rule. In other words, in 
a ID FQ the sequences of the DBRs are fixed by the re¬ 
currence rule. For example, the first five DBR’s are A5, 
BAB^ 5, BAB and BABA. 

In view of the above discussion certain aspects of FQ 
pseudo band gap formation and the nature of the band 
edge states are seen under a new light. In particular, 

• The so-called fundamental pseudo band gaps 
(FBGI and FBG 2 in Fig. p!^f)) are essentially 
parts of a unique band gap separated by the pho¬ 
tonic miniband of coupled cavities, 

• Higher order pseudo band gaps are low- 
transmission regions formed by neighboring 
multi-peaked states of the coupled DBRs (more 
external shaded regions in Fig. [T^f)). The origin 
of their formation is exactly the same as in the case 
of coupled microcavity structures with balanced 
mirrors (see Fig. [T^c)). 

• Fundamental band edge states are of two different 
types. The first are coupled cavity miniband states 
and originate from A A layers. The band edge states 
of second type are essentially coupled DBR states 
and with respect to the coupled cavity states are 
situated on the opposite side of the gap. 

We now examine the properties of the band edge states 
of the first type, i.e. the AA-layers-related coupled cav¬ 
ity miniband of the Fibonacci quasicrystal. As already 
mentioned above, the coupling of cavities inside the ID 


FQ takes place via unbalanced mirrors. This results in 
the formation of photonic states which are not homo¬ 
geneously distributed within the spectral width of the 
miniband, as opposed to the case of a coupled microcav¬ 
ity system with uniform DBRs [2ni[22]. For a more de¬ 
tailed analysis we consider an 8 th order FQ which counts 
eight occurrences of AA-cavities. For simplicity, we re¬ 
name the AA-cavities to Cm (m=I, 2 ,..., 8 ) and rewrite 
the Fibonacci sequence as 

AB Cl BAB C 2 B C 3 BAB C 4 BAB C 5 B Cq BAB 
C7 B Cg B. 


The typical broad range transmission spectrum as well 
as the zoom of the miniband region are plotted respec¬ 
tively in panels (a) and (b) of Fig. 16 As we can see, the 
miniband is composed of three regions of relatively high 
transmission states. While a total of six marked peaks 
can be counted, the number of states involved in the mini¬ 
band is eight whereby two states are less pronounced 
and are not visible in the T{uj) spectrum. The accu¬ 
mulated phase, which is a very useful tool when analyz¬ 
ing unknown resonant structures lailHig], is plotted in 
Fig.jl^c) and counts a total of exactly 87 r when sweeping 
the frequency through the miniband region. The three 
high-T channels are separated by two regions of very low 
transmission. These last are often attributed to the “cat¬ 
egory” of higher order pseudo band gaps, while their ori¬ 
gin is clearly related to the cavity miniband structure. 

In order to understand better the formation of the in- 
homogeneously split miniband we refer to the electric 
field intensity distribution inside the quasicrystal around 
the miniband central frequency ujd (Fig. p^d)). The po¬ 
sitions of the bright spots, which indicate to the large 
filed intensities accumulated in cavities, show certain 
symmetric patterns both through the depth inside the 
structure (x-axis) as well as over the frequency range {y- 
axis) around the uj/ujq = I point. 

As an example, we analyze two of the symmetric pat¬ 
terns, which are shown in panels (e) and (f) of Fig. 16 


First, we notice that cavities Ci and C/i are weakly cou¬ 
pled through the large mirror BABC 2 BCgBAB^ there¬ 
fore the resulting mode splitting is small (on ^-axis). 
More precisely, this double-peaked state is formed by a 
symmetric coupling of three cavities Ci , C 4 and C 7 since 
an identical mirror BABC^BCqBAB is present between 
C4 and C7 : 
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FIG. 16. (a) The broad range transmission spectrum of an 8th 
order FQ. (b) The zoom of the spectrum and (c) the accumu¬ 
lated phase around the miniband region, (d) The scattering 
states map of the miniband states. Panels (e) and (f) ex¬ 
plain the degenerate mode coupling between various cavities 
distributed through the Fibonacci sequence. 


AB Cl [BAB C 2 B C 3 BAB] C 4 [BAB C 5 B Cq BAB] 
Cj B Cs B. 


We notice also that the cavities C 2 and C 3 , which are 
strongly coupled through the very weak mirror 5, there¬ 
fore, the resulting mode splitting pushes the new states 
far from uj/ujd = 1 to the extremes of the miniband 
[Hnni. Moreover, each of these new states itself is a 
doublet state, which can be explained by further analyz¬ 
ing the FQ structure. Namely, as in the case of C 2 -C 3 
coupling, we notice that the same scenario of coupling 
takes place also for C 5 and Cq cavities, which pushes the 
new split states exactly to the same frequencies as for 
the C 2 -C 3 case forming thus a new couple of degenerate 
states at the extremes of the miniband. In this case, how¬ 
ever, the coupling between cavity pairs C 2 -C 3 and C^-Cq 
is much weaker; 


AB Cl BAB C 2 B Cs [BAB C 4 BAB] C 5 B Cq BAB C 7 
BCsB, 


therefore, only a small amount of mode splitting is 
observed. 

Thus, while this kind of analysis can be done further 
considering other cavity pairings (we leave the reader to 
exercise with), the formation of inhomogeneously dis¬ 
tributed miniband states in a ID photonic quasicrystal 
of Fibonacci type is comprehensible now. 


Localisation properties of band edge states 


It is well-known that the energy spectrum of a localized 
at a frequency ujq state is described by a complex function 


t{uj) ^ 


7 

CJ — cjo + iy ’ 


( 12 ) 


and its transmission spectrum has a Lorentzian lineshape 
(Fig.[T7|:a)) 


Here, 27 is the full width at the half-maximum (FWHM) 
of the resonance. In other words, 7 describes how good 
the cavity confines the electromagnetic energy or, equiva¬ 
lently, how lossy it is. Thus, the FWHM in the frequency- 
domain represents the decay rate, proportional to 7 “^, 
of the resonance in the time-domain. In fact, the in¬ 
verse Fourier transformation of a Lorentzian lineshape is 
an exponentially decaying in time function (solid line in 
Fig-lHlb)). 

In case when two identical resonances (same 7 ) are 
degenerately coupled, the resulting lineshape is described 
through a coherent sum of two Lorentzians as 


D{UJ) 


7 

CJ — cjo + 5 + iy 


7 

CJ — cjo — s + iy 


2 


(14) 


where 2 s is the coupling induced splitting of the doublet 
resonance. This lineshape can be also presented as 


D{UJ) 


7 

CJ — cjo — iy 


2 

X 


_ 2 [{U! - U!q)‘^ + 7^] _ 

(cj — cjo)^ — 7 ^ — + 2ij{uj — cjo) 


(15) 


This result has the form L{uj) x in which we 

recognize the single-resonance Lorentzian L{uj) convo¬ 
luted with the function It is easy to see that 

lim D(cj) = I/(cj), i.e. when the coupling is infinitely 

weak, one ends up with a single Lorentzian lineshape. 

In a weak coupling situation (e.g., a strong intracav¬ 
ity mirror) the splitting is small and a clear doublet is 
not formed (dotted line in Fig. [^a)). While this line- 
shape appears as an effectively broader “single-peaked” 
resonance without any evident fine spectral features, its 
Fourier spectrum reveals a clearly faster decay (dotted 
line in Fig. [T^b)) than the pure exponential one of the 
Lorentzian resonance. Moreover, the composite nature 
of this double-state resonance is brought to light nicely 
in the time-domain; it can be observed at a relatively 
longer time scale that the state decay is not monotonic 
but it performs regular oscillations. The time period of 
these oscillations is defined by the inverse of the doublet 
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FIG. 17. (a) Spectral lineshape of a single (solid line) and 
two weakly coupled Lorentzian (dashed line), (b) The corre¬ 
sponding inverse FFT transformations of the single and cou¬ 
pled states, (c) The evolution of the doublet splitting with 
increasing coupling, s/ 7 , (d). The corresponding time decays 
of doublet-states showing the oscillatory-time behavior. The 
insets is the zoom of the time-spectra at the initial part of the 
decays. 

splitting, s~^. This is what typically is observed in time- 
resolved pulse propagation experiments mmmm, 
which essentially reflects the convolution of the ultrashort 
pulse lineshape with that of the composite resonance; 
regular beatings occur between two coupled states which 
manifest as oscillations in the transmitted intensity (see 
Fig.§. In addition, the transmitted pulse is also charac¬ 
terized by a larger delay which is caused by the transient 
time necessary to build up the double resonance. 

Figure [T^c) shows the evolution of the spectral line- 
shape of two coupled Lorentzian states at increasing cou¬ 
pling strength 5/7. The corresponding iFFT spectra, 
plotted in Fig, [l^c), show that the oscillations become 
faster for larger doublet splittings and the intensity de¬ 
cays faster at short times (see the inset of Fig. [T^c)). 

When more than two resonances are degenerately cou¬ 
pled within a photonic structure, the decay of such a state 
may show an even more complicated time-behavior. In 
particular, more than one oscillation period can appear 
in the intensity decay profile, but as well the amplitude of 
beatings can vary because of unevenly distributed degen¬ 
erate resonances within the structure. A photonic struc¬ 
ture which fits well this situation is the ID Fibonacci qua¬ 
sicrystal. As we discussed earlier in this subsection, the 
band edge states of a ID FQ are miniband states originat¬ 


ing from AA half-wavelength layers unevenly distributed 
throughout the photonic quasicrystal. While for their 
localization properties terms such as “exotic” or quasi¬ 
localization are typically used to describe a weaker-than- 
exponential decay, we see that this property is exactly 
what is expected for photonic resonances originating from 
degenerately coupled states. 
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